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The equilibrium behavior of a system of elastic layers under tension in the presence of correlated disorder is studied using 
functional renormalization group techniques. The model exhibits many of the features of the Bose glass phase of type II 
superconductors induced by columnar defects, but may be more directly applicable to charge density waves, incommensurate 
striped magnetic phases, stacked membranes under tension, vicinal crystal surfaces, or superconducting "vortex-chains". Below 
five dimensions, an epsilon expansion for the stable zero temperature fixed point yields the properties of the glassy phase. 
Transverse to the direction of correlation, the randomness induces logarithmic growth of displacements. Displacements are 
strongly localized in the correlation direction. The absence of a response to a weak applied transverse field (transverse Meissner 
effect) is demonstrated analytically. In this simple model, the localized phase is stable to point disorder, in contrast to the 
behavior in the presence of dislocations, in which the converse is believed to be true. 
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Recently, considerable progress has been made in un- 
derstanding the behavior of elastic media in the presence 
of randomness. Examples include single flux lines in a 
dirty superconductor , interfaces in random magnets 
[^-|| , charge density waves J0| , and the vortex glass phase 
of bulk superconductors The experimental work of 

Civale et. al. ^ has demonstrated the feasibility of creat- 
ing superconducting samples with correlated (columnar) 
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FIG. 1. A stack of layers in three dimensions fluctuating 
in the presence of correlated disorder. Dark and light thick 
vertical lines indicate individual correlated pins, in the fore- 
ground and hidden from view, respectively. The spacing along 
the layering (a;i) axis, I, and the displacement fleld u{x,z), 
are indicated. 

disorder. One of the most striking aspects of the result- 
ing localized Bose glass phase is the existence of 



a finite critical mismatch angle 6c between the applied 
field and the correlation direction, such that the fiux 
lines remain parallel to the correlation axis for 6 < 9c- 
In this paper, an analogous localized phase is studied 
analytically near 5 dimensions for a somewhat simplified 
model which may describe other systems of interest. The 
puUing-away from the correlation axis for 6 > 9c is dis- 
cussed in Ref . Q . 

Consider a system made up of a stack of layers (for 
general d, these will be oriented manifolds of dimension 
d — 1, for instance interfaces) which may fluctuate in the 
perpendicular direction, but may not pass through one 
another (see Fig.Q). Such a model could describe a charge 
density wave, the domain walls of an incommensurate 
striped phase in a magnet ||l3| , a stack of membranes un- 
der tension, or, for c? = 2, a set of steps on a miscut crys- 



tal surface 1 14|. Another possibile realization is suggested 
by recent observations of "vortex-chains" in YBCO , 
which should fluctuate as layers defined by the chain and 
magnetic field directions. The displacement of the A:*'' 
layer, Uk, is defined by Xk,i = k£ + Uk{x2, ■ ■ ■ , Xd-i, z), 
where £ is the average layer separation, and xi is taken to 
be the layering axis. The last coordinate Xd = z has been 
distinguished as the direction of correlations for the ran- 
dom potential. The use of a displacement field neglects 
dislocations, which will be important in many systems 
beyond some length scale. Taking the continuum limit 
Uk{x2, . . . , Xd^i, z) u(x, z), the Hamiltonian is 



H = Jd'^-hcdz\^^ (V 



hdzU 



+Vc{u,^) + Vp{u,^,z)}, (1) 

where the xi coordinate has been rescaled to remove the 
anisotropy of the in-layer elasticity, and a momentum 
cut-off A ~ l/£ is implicitly included in the x direction 
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(but not in z) due to the discreteness of the layers. A 
non-zero h represents a force tending to tilt the layers, 
such as that caused by a change in the applied field in a 
superconductor, or by tilted boundary conditions. Only 
the response to a small h will be considered here. Vc 
and Vp are random potentials describing, respectively, 
columnar and point disorder. By including both types of 
pinning, it is possible to study the competition between 
vortex-glass and Bose-glass like phases in this dislocation- 
free model . Because of the periodicity of the stack 
of layers, these potentials must be periodic functions of 
u. For weak disorder, the fixed point potentials will have 
Gaussian distributions, with the correlation functions 

{Vc{u,^)Vc{u',^')) = Rc{u-u')Si^-^'), (2) 
(yp(M,x,z) Vp(ii',x',z')) — Rp{u — u')S(x~x')S{z — z'). (3) 

Choosing £ — 27r, periodicity implies Rc,p{u + 27r) = 
Rc,piu)- 

To understand the behavior on long wavelengths, a 
renormalization group (RG) analysis was performed us- 
ing the methods of Ref. Q . Lengths are rescaled parallel 
and perpendicular to the disorder, according to x — > 6x 
and z — > b^z. The displacement cannot be rescaled due 
to periodicity. The cut-off is kept fixed by integrating 
out modes with x momenta in a shell A/6 < p < A (mo- 
menta in the z direction are unrestricted). Formally, this 
procedure can be carried out by expanding the partition 
function in Vc and Vp and performing the functional in- 
tegrals order by order. The resulting terms arc functions 
of the remaining modes, and upon reexponentiation yield 
flows of the interactions in the Hamiltonian. Since u is 
dimensionless, however, Rc and Rp cannot be taken to 
have the form commonly encountered in RG studies 
of phase transitions. 

For d = 5 — e, the behavior of the system is dominated 
by a non-trivial zero-temperature fixed point. The scale 
changes yield the eigenvalues At = —d — —2 — X + 0(e) 
and Aflp = — x -|- 0(e). Naively, x = 1 + 0{e), but this 
point will be returned to later. Regardless of the value 
of X (so long as it is positive), T and Rp are formally 
irrelevant, so I will begin by working directly at T = 
Rp = 0. 

At zero temperature, the computation of the partition 
function reduces to the optimization of the hamiltonian. 
At each step of the RG, the lowest energy configuration 
of the modes in the shell is found as a function of the low 
momentum modes, which are held fixed (see, e.g. Ref. 
Q). The minimum of H (Eq.|l|) clearly satisfies dzU — 
exactly. The renormalization of Rc, therefore, must be 
identical to the case of point disorder in d — 1 dimen- 
sions. In addition, the statistical Galilean invariance of 
the d—l dimensional model [ pT[ guarantees that temper- 
ature is only trivially renormalized by the scale changes. 
There is, however, a non-trivial renormalization of K, 
since the iterative minimization does not throw out ex- 



citation information until it is on a scale smaller than 
the cut-off. Defining the force-force correlation function 
Ac{u) = —Rc{u), the RG equations for b — e'^^ to lowest 
non-trivial order are 



dK 

— =(2-2x-Ae(0)/(8^^ 



dAc 
IT 



1 

8^ 



{A'cf + ^'c (Ac-Ac(O)) 



(4) 
(5) 



where, for simplicity, I have taken K = 1. Eqj^ is ob- 
tained assuming analyticity of Ac, and will be corrected 
later. Eq.|^ has been derived previously by a number of 
other authors p^,|6|,p9t , and has the stable 27r-periodic 
fixed point solution 

Ac(") - [{u - - ttVS] , for < u < 2^ . 

(6) 

This fixed point solution leads [|l9| to logarithmic dis- 
placement fluctuations. First order perturbation theory 
in Rc, evaluated at the fixed point (Eq.||), gives the 0(e) 
result 
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(7) 



where u(x) is defined as the z average of u(x, z), and the 
angular brackets denote both thermal and disorder aver- 
aging. The fixed point function has a slope discontinuity 
when M = 0, so that A(^(0) = — oo! Note that the proce- 
dure employed in the near-threshold dynamic problem of 
Ref. [|j for handling this singularity does not apply here 
due to the different physics of force-free equilibrium. The 
divergence implies that the feedback of Ac {u) to the elas- 
tic terms (Eq.^) must be re-analyzed. With no assump- 
tions on the analyticity of Ac(u), the term generated by 
the RG which led to this equation takes the form 



dH 
IT 



A3 



clastic 



1 

2l67r2i^i/2 
A 



X exp 



^1/2' 



d'^-hidzdz' /:\c{u{z)-u{z')) 
(8) 



The exponential decay for large separations justifies a 
gradient expansion of u{z) near z', which combined with 
the smah u behavior Ac(m) w ^{n'^ /i - 7r|u| + u'^/2) 
yields a new term in the Hamiltonian 



Ai? = I Jd'^-hcdz\d,u\. 
Eq.^ is replaced by the pair of equations 



da 
H 

dk 
la 



(2-x)'t + 



7rAeXi/2 



(2-2x-e/3))i^, 



(9) 

(10) 
(11) 
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The new elastic term has a direct physical meaning. 
Consider the response of the system to an infinitesimal 
field h. In the absence of disorder, the minimum of the 
Hamiltonian as a function of dzU is shifted over by an 
amount linear in h, resulting in a response 9 = dzU oc h. 
This will be true regardless of the magnitude of K. If the 
elastic term has the form of cr|92u|, however, the mini- 
mum will remain at dzU — for all h < K. This implies 
the existence of a finite threshold force, below which the 
layers remain locked in their localized positions. 

It is clear, both from the physical interpretation above 
and from the form of Eq.|l^ and Eqjl^, that the sim- 
ple expectation of 0(e) corrections to x "^iU not hold. 
Eq.|l^ suggests a value of x = 2, in accord with a simple 
"random walk" picture of excursions from the localized 
ground state. A careful analysis, however, must continue 
the RG procedure after the new elastic term has been 
generated. While this cannot affect the renormalization 
of Ac, it will probably affect the z-dependent portions 
of the Hamiltonian. It is not known whether a consistent 
treatment of such a non-analytic elastic term is possible 
even to 0(e). Justification of the value % = 2, even to 
lowest order, requires a detailed investigation of the zero 
temperature RG along the lines of that performed in ap- 
pendix C of Ref. . Such an analysis is not available at 
present. 

Because the vanishing of the tilt response requires 
Ap(0) = — oo, the weak field behavior is sensitive to cor- 
rections from formally irrelevant operators such as tem- 
perature and point disorder. To safely conclude the sta- 
bility of the correlated phase, the corrections will now be 
studied in more detail. These operators yield additional 
renormalizations of Ac 



dAc 



dl 



A 



T,Af 



167r2i^i/2 



A^T 



1 



Ap(0) 



(12) 



Because of the presence of such terms, analytic behav- 
ior of Ac{u) will persist within a narrow boundary layer 
around u = 0. Since temperature and uncorrelated disor- 
der are irrelevant operators, the width w of the boundary 
layer decreases under the RG (note that since 6 > \Xrp\, 
point disorder will dominate this width at long length 
scales). Although this rounding is not a property of the 
fixed point value of Ac, the rapid divergence of Aq{u) 
in the absence of the irrelevant operators (see, e.g. Ref. 
1^ ) indicates that the size of the smoothed region is de- 
termined by the terms in EqJT^. 

As the RG iterates to longer length scales, Ac{u) 
sharpens up near the origin. From Eq.^, the jump 
in slope across the boundary layer, where Ac{u) must 
match its fixed point value, is 0(e). Since this change in 
slope must be accommodated over a width w, the cur- 
vature Ap(0) —e/w. Equating the terms in Eq.[l^ 
to those in Eq.|| gives the scaling of the boundary layer 
width 



;(Z) 



A 



A'T{l) + -Ap{Q;l) 



(13) 




FIG. 2. A schematic illustration of the columnar disorder 
correlation function Ac{u). For \u~2nk\ 3> w, Ac{u) has the 
simple form of Eq.6. For \u — 2Tvk\ <C w, Ac{u) is rounded, 
with curvature of 0{e/w). 

Ac{u) is illustrated in Figj|. From the structure of 
Eq.^, it is clear that (for small gradients) the generated 
terms take the form 



dH 

IT 



clastic 



eA^ / d'^-^yidz f 



A 



\dzu\ 



(14) 



where f{u) \u\ for w <C |u| <g; 1, but is smooth for 
|u| ^ w. 

Tilting for small h can only occur for 9 within the 
smooth boundary layer. Although Eq.|l^ is valid only 
while feedback from a can be neglected, it does demon- 
strate that this boundary layer width, w, vanishes ex- 
ponentially as a function of length scale. Since the av- 
erage tilt is a bulk (q = 0) quantity, the boundary 
layer does not contribute. By contrast, it is precisely 
the small boundary-layer in Eq.|lJ which leads to ther- 
mal creep in weakly driven interfaces and charge-density 
waves |20(] . The only important effect of the interaction 
between point and columnar disorder in this model is to 
decrease the relevance of point disorder, both through 
an increase in x and through terms from Ac feeding into 
the RG equation for Ap (not shown here). 

In conclusion, the system of elastic layers studied here 
forms a disordered-dominated phase in the presence of 
correlated randomness. Despite the reduction in the 
number of components of the displacement field and the 
neglect of dislocations, the model exhibits analytically 
many of the properties of the Bose glass. The physi- 
cal properties are governed by a zero temperature fixed 
point, at which the layers are completely parallel to the 
direction of correlation. Due to this localization, the tilt 
response to an applied transverse field vanishes below 
some finite threshold field, corresponding to the "trans- 
verse Meissner effect" of Ref. The fluctuations of 
the layers perpendicular to the correlations grow loga- 
rithmically, with a universal prefactor. In contrast to the 
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Bose glass case the system is stable against uncor- (1992). 
related disorder. It thus appears that dislocations are 
necessary to make point disorder relevant in the local- 
ized phase. Naive arguments suggest a simple "random 
walk" scaling, x = 2, for the low-lying excitations, al- 
though it is quite likely that this result will be corrected 
by non-trivial renormalizations. The analysis of these ex- 
cursions from the ground state and the extension to more 
components remain interesting open problems. 
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